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1 Introduction 



The realm of integrable systems in 2d statistical physics or (l+l)d QFT has attracted much 
attention for a long time. Important contributions in this development were the solution of 
the anisotropic sl{2) Heisenberg model (XYZ-model) by means of the Bethe ansatz [Q and 
its subsequent reformulation in an algebraic version, the algebraic Bethe ansatz or so called 
Quantum Inverse Scattering Mathod (QISM) 0. Generalizations to higher rank groups, which 
allow to treat models with internal degrees of freedom, were in this context achieved in . The 
efforts resulted in a plethora of models soluble by means of the QISM, insofar that eigenvalues 
and eigenvectors of the Hamiltonian were found explicitely. This method also stimulated the 
investigation of various areas in mathematical physics, such as quantum groups, theory of knots 
etc. Despite the indisputable achievements of the QISM and the rather simple action of the 
inverse problem operators, which can be interpreted as creation and annihilation operators for 
quasiparticles and as generating functions for the conserved quantities respectively, the study 
of correlation functions and formfactors has proven to be rather intricate. This is partly due to 
the fact that the solution of the inverse problem (expressing the original microscopic operators 
by means of the operators figuring in the algebraic Bethe ansatz) has only been achieved 
recently [^, |^, 0, and by the fact that the action of the quasiparticle creation and annihilation 
operators which figure in the construction of the eigenfunctions (Bethe wavevectors) is obscured 
on the level of the microscopic variables (spin raising and lowering operators) by nonlocal 
effects arising from polarization clouds (compensating pairs of local raising and lowering spin 
operators). In a seminal paper Maillet and Sanchez de Santos revealed an application of so 
called Drinfel'd twists to obtain a basis for the sl{2) XXX and XX Z model which allows to 
express the creation and annihilation operators in a completely symmetric way with the further 
advantage of being polarization free, that is being built from the respective quasiclassical 
Gaudin operators dressed diagonally, thus supressing non-local effects for these operators. A 
generalized transformation has subsequently been applied to the sl{2) XYZ model and 
been used to resolve the nested hierarchy in the Bethevectors of the sl{n) XXX model |p. 
In this paper we report the construction of suitable Drinfel'd twists for the Z„ x Z„ symmetric 



Belavin model |T0[, which can be thought of as a n-state generalization of the sl{2) XYZ 
model. Our results generalize the findings for the sl{n) XXX model and provide a completely 
symmetric representation of the creation operators as well as a resolution of the hierarchical 
structure of the Bethe wavevectors. 

The paper is organized as follows: Section 2 provides a short survey of the Belavin model and 
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its reformulation as an face-type model. Section 3 deals with the factorizing twists and the 
computation of the operator valued entries of the monodromy matrix. In section 4 we discuss 
the Bethe wavevectors and section 5 contains the conclusions. 



2 Belavin model and the corresponding IRF model 



A possible n-state generalization of the eight- vertex model |]T| is given by the Z„ x Z„ symmetric 
model of Belavin , whose Boltzmann weights fulfill the Yang-Baxter equation (YBE) 



S{z, - z^^S. - ^3)S£ S{z, - Z3)JJ = S{z, - z^^^^l^ S{z, - z^^;^ S{z, - J (1) 

(here and subsequently double indices k mean summation over 0, 1, ... ,n — 1 unless stated 
otherwise). It can be parametrized in the following way [|1^] [| 



with 



n—l /n— 1 



h{z) = l[e^^\z)[l[e^^\o)] (3) 



1=0 \i=l 



where 9^^\z) = represents the theta function of rational characteristics 1/2 — 

2 

i/n, 1/2. The theta function of characteristics a, 6 G M is given as a Fourier series (r is a 
fixed complex number in the upper half plane, and A,- = Z + Zr is the lattice generated by 1 
and t) 



^ ^ meZ 



) (2;) = exp (7rir(m + a)z + 27ri(m + a){z + b)) (4) 



and has zeros at z = (| — 6) + (| — a)r mod A,-. 

The matrix (0) is unitary, S{z)S{—z) = I and obeys the initial condition S'(O) = P [|1^, where 
P is the permutation operator. 



There exists a vertex- face map which transforms the Belavin model into a face type model [14 



which in turn can be thought of as a multicomponent generalization of the six-vertex model. 
^The slight difference compared to originates from the normahzation of the S-matrix. 
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The correspondence is given by 

Siz, - z^)Mli{z,) ® M:^^{z^) = M„^;,,(zi) ® M:'(z2)i?(^i - ^2|a)p> (5) 



where R{a\z) has the form 



h{z + w) ^s" - + w{a - z/)'^'^ + w) fj^ giu{a + fi) fj^ gi^{a) 

f^,u( I N _ + ^) H^'^ - + w(a - i>Y''' + w + z) 
^ ~ h{z) his" - s^' + w{a - uY''^ + w) 

( 1 if z > /i 

g^a) = I - + wa°''^) ifi = (8) 

[ — + wa^'^ — w) if < i < ;U 

where are arbitrary complex numbers, a G ^^Zq CA^ with weights of the affine Lie 
algebra A^n-ii A = ~ = 0, 1, . . . , n — 1 and Aq = A„. The a^'^ is given in fl^ and 

obeys (a — z>)^'^ = 0. 

The intertwining vector Ml^{z) is given by 



n-l 



= *(0:^Wo,...,0^:Wn-i)n^?^(«) (9) 

1=0 

(j)'^{z)j = e^^\z-nwa^'''^). (10) 



Applying a sum of threefold tensorpro ducts of functions (|) to (P using the vertex face map 
and exploiting the explicit index structure of (||) we obtain a modified Yang-Baxter equation 

R{\z, ~z,\a + a,);\2R{z^ - z,\aY^\'-;R{z, -z,\a + ^,^1^ 
= Riz2 - z,\aYl:';^|Riz^ -z,\a + fi2)Z2R(^^ ' ■ (^l) 

This can be written symbolically as (we set R{zi — Zj\a) = Rij{a)) 

Ruia + (5-3)i?i3(a)i?23(a + ^i) = i?23(a)^i3(« + '3"2)^i2(a) • 

The monodromy matrix for a chain with sites associated with the R-matrix (^ is 

7'o,i...Ar(A|a) = i?o7v(a + ai + . . . + ctn-i) ■ ■ ■ Ro2ia + ai)Roi{a) (12) 

where denotes the horizontal auxiliary space with spectral parameter A; 1, . . . , A^ label the 
vertical quantum spaces (with local inhomogeneities {zi}) whose tensorproduct constitutes the 
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physical Hilbertspace TiN- An algebraic construction of eigenvalues of the monodromy matrix 
associated to the matrix (|2D was performed in [0 and parallels the procedure for the eight- 
vertex model in (One has to take into account that our construction of the monodromy 
matrix (|12D differs from the monodromy matrix in |jl6[ by an additional transformation in the 
physical space). 

Subsequently we will focus on the construction of a factorizing F matrix for the R-matrix 
and pursue its consequences for the structure of the monodromy matrix (|12|). 



3 The F basis 



The factorizing F- matrix for N sites (A^ quantum spaces), being defined by the relation 
Ri...n{o) = -^(7(1 Ar)('^)-^i...Af('2) 0, turns out to be given by formally the same expression 
as found in for the sl{n) XXX- model 

* N 

where P"' projects on the ttj-th component in the i-th. space and the labels ao-{j) satisfy the 
conditions 

oia{i+i) > a<7(i) if + 1) > (t{{) 

> a^(i) if cr{i + 1) < ^(i) . (14) 

The modification of the Yang-Baxter equation ([11]) requires a particular rule for the handling of 
the parameter a in the formation of the intertwining matrix R"{a) (related to the permutation 
0"), which can be read off from the modified composition law 

R^'^^a) = i?.(,),.(,+i)(fi,)i?'^(a) 

Oi = a + + . . . + 6-^(i_i) (15) 

where cTj is the transposition oii,i + 1, and a an arbitrary permutation. 
R^{a) has the intertwining property 

i?'^(a)To,i...iv(a) = To,<^{i)...a(7V)(a)P''(a + ^h) (16) 

where ah is associated with the matrix indices in the space 0. 
The matrix Fi jv(a) satisfies the factorizing equation 

RL.M = F;^l ^^{a)F,,„r,{a) . (17) 
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A proof of the latter equation can be found in and relies on the fact that 



PtP"R^J = PtP^'hj . (18) 

The modification of the composition law induced by the parameter a being insignificant. 
We will proceed by computing elements of the monodromy matrix in the new basis provided by 
the F-matrix. The transformation law To,i...Af('2) = Fi...Nio.)To^i,,,N{o.)F^^j^(a + ah) is enforced 
by the requirement that the resulting operator is sjTiimetric, i.e. To,i...Ar(a) = To_o-(i)...o-(Af)(fl)5 
which then follows from (|16|). 



The computation of the diagonal element T„_i„_i(A|a) proceeds along the same lines as that 
in the generalized sl{n) XXX-model in 0. Let us consider the action of the matrix F{a) on 
Tn-in-iio,): 

* N 

Fi...iv(a)T„_in_i(a) = E E n^"«'^i...^(«)^o""'^o.i-^(«)^o'^~' 

* N 

= E E n^r(?^o"'%^(i)...<x(^)(«)^o"~'^i...7v(a + ^/.)(19) 

The specialization to the entry {n — l,n — 1) of the auxiliary space here is achieved by the 
projectors Pq ~^- For the second equality in ( [TP| ) we have used relation (|TB|) and the obvious 
fact that Pq~^ commutes with j^f{a)- To simplify the following argument we distinguish in 
the sum Yl* cases of various multiplicities of the occurrence of the group index n — 1: 



N N N-k 

' ^(i) 



o-eS'jv k=0 a^(^iy..a^fff) j=N-k+l j=l 

xPo"-iro,.(i)....(^)(a)P^-ii?^..,^(a + an) . (20) 



Let us consider the prefactor of R1 j^{a + ah) on the r.h.s. of Eq. ( pOD more closely. Using 
specific features of the i?-matrices we can rewrite it as follows: 

N-k N 

WCi^' n P:u]Po-'T,,.W.MN){a)Pr' 

j=l j=7V-fc+l 

N-k / N-l X n-ln-1 

= n ^"o?' I ^o,a(Jv)(a + '^'^m) 1 • • • {Po,a(N-k+i){a))l_^l_^ 

i=l V 1^=1 / n-ln-l 

N 

xPr'To,^ii).MN-k){a)Pr' n P^U) 

j=N-k+l 



N-k N 

nC)'^o"-^To,.a)....(^_.)(a)Po"-^ n 

j=l j=N~k+l 

N-k / i-1 \ "-l'"<T(i) N-k N 

n ( Ro.i^) {a+Yi ^^im) ) ) n p^s' n p^ui 

j=N-k+l 



(21) 



m=l 



Inserting the r.h.s. of (|2TD into Eq. ( pOD one sees that the product Yli yRoa{i) (ct + Ylm=i ^a-(m)) 

creates a diagonal dressing factor for T„_i„_i(a) and the product of projectors apphed to 
R'^(a + ah) gives Fi j^{a + ah)- One obtains 

f„_i„_i(A|a) = Fn-i(a) diag{6(A - Zi), . . . , 6(A - z^), 1} . 



-l."<T(i) 



(22) 



with the abbreviations 

j-l ( n-l 



n n 

m=0 L j=0 



Qij iflm) 



n 



, . 

-TT gimjam + (jm - 1)^) 

7=0 fi'im(a^ + (jm - l)rh + j) 



X 



/i(s™ - + ^[S;; + (j„ - 2)m]-'^') 



6(A) 



/i(s'" - s^' + w[am + {jm - 2)m]™J + w) 

m— 1 

a + ^ ki i 

1=0 

h{X) 



h{X + w) 

and km gives the multiphcity of the upper matrix labels = m. 



(23) 



To compute T„_i„„2(A|a) one has to distinguish in the sum J2* cases of various multiplic- 
ities kn-i and kn-2 of the occurrence of group indices n — 1 and n — 2: 

N N-kn-i *" AT N-k„-i 

F,...N{a)T.-.n-.{a) = E E E E 11 ^."(^0 11 ^^ilL) 

O-GSjv fc„_l=0 fc„_2=0 ao-(l)---"o-(JV) in-l=A^-fcn-l + l in-2 = A'-fcn- 1 - fe„-2 + 1 

N-kn-l-k„-2 

X n ^ra^n""' V(i).-w(«)^o""'^L..iv(« + ^h) . (24) 

Evaluating the matrix product in To,o-(i)...cr(Ar)(a) leads to 

AT N-kn^l N-kn-l-kn-2 

n ^"(i.) n n c,^fprTo,.i^).MN){a)pr 

J„_l=7V-fc„_i+l j„_2=iV-fc„_l-fc„^2 + l i=i 

n-ln-2 

E (^oiv^i(^^-fc-i))„_i„_2--- (^om(«^+i))IlIr2 

i=7V-A;,i_i-A;,i_2+l 
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k=l jj^i 
N-kn-i N-k„-i 

j=A''— fe,i_i — fe,i_2+l l=i+l 

N-k„-l-kn-2 

X n &"~''°KA-^^iafc)i?u„_in^r^""'^"-i,-2 (25) 

fc=l jVi 



with the abbreviations 



m=l 

i = a(i) 



One notes that in the calculation the index a-^ has changed from n — 2 to n — 1. As the 
distribution of a's is therefore no longer consistent with the conditions (14) in the sum ^* 



one has to correct this fact by commuting the site i through all higher sites j with a-j = n — 2. 
So, taking into account (0), one has to insert an additional factor 



R-^n-" {O'j.n-^ + o'h) ■ ■ ■ R-'^~rr{or^+ <^h) (26) 

iN-kn^l^ N-kn-1 ' ll+l^ l + l ' \ ' 

between the projectors and ^{a + au) in (|^). Because of Eq. (|l^) no further corrections 
are neccessary. For the following calculation two equations are needed 



and 

fpi pn— 1 p Tpj r>n-l fnc 

with {Pi = Sai,/3jSaj,i3i the permutation operator in space i and j. 



Let us now concentrate on the term with i = N — kn-i — kn-2 + 1 in (|25|) and use 
to create the needed factor (p6l). Because of (|^) the second term in (^Tl) gives rise to an 
El^_2n-i with i 7^ i. So the only possibility to get an -E*_2„_i is to use the first term in (^). 
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Corrections in the other terms with j > N — — fc„_2 + 1 cannot lead to an expression with 
-En-2,n-i j has not to be commuted with the site i. So the only term that contains -E^_2_„_i 
after the corrections for that special R1 ^ in (|25|) is 



N—k„-i ,n—ln—2t\ 1-^ \ N—k„~i—k„-2 



1=1+1 ^ * " fc=l 

N-kn-i ^ N~k„-i~k„-2 

= n fof^^l^^""''""'^^-^^!^"-^) n &""''"HA-^da.)EL2,n-i- (29) 

i=j+l ^ * ''^ A;=l 

Because of the symmetry of To^i...Ar(A|a) all other terms have to be of the same form as (p9|). 
After taking into account the projectors the resulting expression is 

N 



T„_i„_2(A|a) = F„_2(a)^c"-2'"~i(A-^,|a„_2)K-2,n-i 

i=l 

®f^, diag{6(A - zj), ...,biX- zj), b{X - z,)b-\zi - zj), 1}. (30) 



For the calculation of T„_i„_3(A|a) one has to distinguish the cases n — 1, n — 2 and n — 3 
in the sum ^2* . The only difference compared to T„_i„_2(A|a) is a term containing a product 
El^^^,n-2^n~2n-i ^^^w showiug up in the matrix product in To^o-(i)...o-iv(a)- One again has to 
correct the distribution of a's with the analog of the equations ( p7| ) and ( ]28|) 









a) 


l)Qi,Qj 


[Zi 


-^3 


a) 



El,PfP,, = Ei,Pt (32) 

and also with a new relation which has to be taken into account when dealing with the term 
containing the product ^;_3,„_2 ^^2,™-! : 

En-3,n-lE'j -Pj.j = -E„_3^„_2-£'n-2,n-l- (33) 

This reasoning leads to {bik = b{zi — Zk)): 



N 

f„_i„_3(A|a) = F„_3(a)^c"-3'"-i(A-^,|5„_3)^;-3,n-i 



i=l 

®f^idiag{6(A - Zk), . . . ,b{X - z^), b{X - Zk)b7j^, b{X - Zk)br^, 1} 

N 

+F„_3(a) /2(^ - 3, n - 2; n - 1)£;;_3,„_2 Ei_,^^^, 

®£,,,^.diag{6(A -Zk),...,b{X- Zk), b{X - Zk)b-l b{X - Zk)b-l, 1} (34) 



with 

/2(n-3,n-2;n-l) = c"-2'"-i(A - z,|5„_2)c"-''"-'(A - z.|5„_3) 

[6«-i'"-2(A - z,|5„„2)c"-''"-^(A - Z,|fi„_3) 

(35) 









an-2) 


5n 






an-2) 



where the second term in (l35| ) has its origin in the second term in (0) and in ( p3D 
Proceeding in an analogous manner one obtains in the general case 

a 

fn-laiMa) = E E E 

k=l n— 1— a<ni<...nfc<n— 1 

^£1. ® E^£n, ® • • • ® ^tLi /'^(^i' n2,...,n,;n- 1) F„_i_«(a) 



diag \ b{\ - 2;^), . . . , 6(A - 2:^)6^ J, . . . , 6(A - . . . , 6(A - ^jO&iJ-, . . . , 6(A -Zj)b^, 1 

(36) 



"2— «i n— 1— rife 



with /fe(ni, 77-2, . . . , TT-fc; tt-o) being defined recursively by 



/fc(rai, ?i2, • • • , n^; no) = - , noTn^ ^0%"' /fe- 1 (^1 ' ^2 , • • • , "-o) + Co^'f Vfc-i('^i, ^^2, • • • , nfc-i; nfc)(37) 



where 

/i(ni;no) = Co'4"° 



Remark : As the above procedure to obtain the operators of the monodromy matrix in the F- 
basis did not rely on the invariance of the monodromy matrix under combined sl{n) rotations 
in the auxilliary and quantum space, the sl{n) XXZ- model can be treated in the same way, 
and we obtain expressions for the T^-i a similar to that in (p6|), with the diffenence that the a 



dependence vanishes, and the parametrization of the elements of the R-matrix is trigonometric 
instead of elliptic (c(A) = J^tx+l) ^ ^(^) = ./nMA+l) )- 



^In the rational sl{n) case this recursion relation leads directly to the result in Eq. (42) in ||]. 
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4 Bethe wavevectors 



Our presentation of the hierarchical Bethe ansatz will be rather sketchy. For details we refer 
the reader to |T^|. 

The operators Tn-ia {a < n — 1) act as quasiparticle creation operators and satisfy the 
Faddeev-Zamolodchikov algebra 

Tn-lMa)Tn-lp{^l\a + a) = T„_i^(/i|a)T„_i5(A|a + ^)R{\ - /i|a)^]j (38) 

or in components of the R-matrix 

[T„_ia(A|a), T„_ia(/i|a + d)] = (39) 
T„_iQ(A|a)T„_i/3(/i|a + d) = ^^^^^ _ ^^^^ T„_i/;(/i|a)T„_ia(A|a + jS) 

c^"(/i- A|a), 



6/3°(/i- A|a) 



T„_i/3(A|a)T„_i,(/i|a + /?). (40) 



Inspired by [^], an ansatz for a Bethe vector is given in terms of a linear superposition of 
products of operators Tn-ia acting on a reference state l^j^'': 



n-l 

'AT 



^„(Ar;Ai,...,Ap|a) = J] ^n-lai(-^l|fl)2Ti-la2('^2|<^ + di) . . 

(41) 



where the reference state is constituted as a A^-fold tensor product of lowest weight states 



.;i^)=*(0,...,0,l) in Ci^) 

and the ^ax,...,ap denote some c-number coefficients, which themselves have to be chosen s.t. 
they are components of a sl{n — 1) wavevector, leading to a nested structure finally giving a 
s/(2) eigenvalue problem. 

It is important to note that the reference state is invariant under the F-transformation 

N — A'" 

due to the special form of the R-matrix. 

We will not impose the Bethe ansatz equations for the spectral parameters {Aj} which turns the 
vector (|4l|) into an eigenvector of the transfer matrix, that is we consider the Bethe wavevector 



being "off-shell" |17 
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In what follows we want to determine the functional form of such vectors, using the explicit 
form of the operators relevant for the Bethe wavevectors. 

In order to clarify the arguments employed in the course of the computation we will present 
the case of s/(2) and s/(3) in quite a detail. The generalization to the general case of s/(n), 
n > 3 will then be rather straightforward. 
For the s/(2) case we have from ( pOD 

Tio(A|a) = f;c°'^(A-.,|a)a«®,^,(^(\-"^-)^^' (42) 
i=i ^ ^ b1 

^2{N;Xi,...,Xp\a) = Tio(Ai|a)Tio(A2|a + 0) . . . Tio(Ap|a + (p - 1) • 0) ^SJ^ 



= J2 i?fnAi,...,A,|z.„...,^.Ja)aJ^)...aJ'')fig). (43) 

The c-number coefficient B^'^\{\i}\{zi}\a) is, due to the "diagonally dressed" spin raising 
operators o"!^ in (^2[), of the form 



B^^\\i, . . . , Xp\zu . . . , Zp) = ^ fJc°'i(A„-2;,(„)|a+(m-l)-0) 



b{Xm -~ Za 



(0. 



(44) 

We now turn to the sl{3) case. The strategy will rely on the symmetry of the wavevector under 
the exchange of arbitrary spectral parameters (the verification of this fact follows the same 



lines as the one in |jT8| using (|TT]) and (|38|) ) which enables us to concentrate on a particularily 
simple term in the sum (PD, and the repeated use of the Faddeev-Zamolodchikov algebra. 
These ideas lead us to propose the following form for the Bethe wavevector 

^3(^1 Ai, . . . , Xpg] Apo+i, . . . , Apo+pja) 

= ^ B^JiKo + l^ ■ ■ ■ ^ Ko+Pi\^(t(1), ■ ■ ■ , Xa{p^)\a) 
Pi po 

X J] J] b'^'{\,^k)-\ail)\a + {k-l)■6 + {l-p^-l)■i)-' 
k=i i=pi+i 

X ^2i(Aa(pi+i)|a) • ■ ■T2i{Xa{po)\a + {po -pi - 1) ■ i) 

X T'2o(A^(i) |a + {po ~pi) ■!)... r2o(A^(p,)|a + {po - Pi) ■ 1 + {pi - 1) ■ 0) n^^^ . (45) 

Consider a special term in the sum (^) of the form (which is motivated by the fact that the 
associated coefficient $ is especially simple to compute, see below) 

T2o(Ai|a)T2o(Ai|a + 0) . . .T2o(ApJa + {pi - 1) ■ 0) 
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X f2i{Xp,+i\a + Pi ■ 0)r2i(Ap,+2|a + Pi ■ + 1) . . . T2i(ApJa + pi ■ 6 + (po - Pi - 1) ■ 1)^5?^ . 

(46) 

Commuting all T2o(Ai|a) to the right using the first term in (|40|) yields an additional factor 

pi po 

n n {&''°(A, - A.|a+(x-l)- 6 + (y-pi-l) 

x=l y=pi + l 

It has to be noted that the associated c-number coefficient $a.^.oi...i in ( ^61) is not evaluated in 
the s/(3) F basis. It can however be expressed in the form (^11) as it is invariant under the 
action of the sl{2) F- matrix. This is due to the fact that it constitutes a component of the 
sl{2) vector whose labels (a non-decreasing series of with respect to the original ordering of 



sites i) correspond via (jM]) to the identity permutation in the definition of the F- matrix (13). 

Invoking the exchange symmetry we arrive thus at the formula (|i5D . 

The creation operators with respect to the lowest weight state are of the form (cf. (|36|)) 

N 



r2i(A|a) = yi(a) J2 c'''(A - ^i|ai)E« diag{6(A - z,), 6(A - z,)6-.\ l}y] (47) 
1=1 

TV 

T2o(A|a) = ^ c°'^(A - 2;i|ao)Ej^ O^^i diag{6(A - Zj)bij^, b{X - Zj)b^^, 
1=1 

+ terms involving e[^1 Elj'l . (48) 

The second term in ( ^81) annihilates the vacuum fij^^ , which is why we did not cite the explicit 
form of the prefactors accompanying these roots. The form of the creation operators permits 
us to further simplify the wavevector (^5]). Taking into account the action of the roots on the 
respective dressing as well as the fact that both groups of creation operators generate a factor 
similar to the s/(2) problem, we obtain 

^3(^5 Ai, . . . , Ap(,; Apo+i, . . . , Apo+pj \a) 

Pi Po 

= E E n n ^'•°(A.(/)-A.(,)|a + (fc-l)-0 + (/-pi-l)-i)-U(A.(0-^.J 

ueSpQ ii^...j^ip^ k=l l=pi+l 

^ -Spo-pi('^'^(pi+i)' • • • ' A(T(po)kjpi+i) • • • ) ^ipo\^) -^pl^(Apo+i7 ■ ■ ■ 7 Ap(,+pJA^(i), . . . , Xa{pi)\a) 
X Sjf (A.(i), . . . , A.(p,)|z,,, . . . , z,^^ \a + (po - Pi) ■ i) i^ij^^^ . . . E^^^E^;^ . . . i^Jj^fig^ 

(49) 

where we defined generalization of (HI) 
-Bp'^(Ai, . . . , Ap|2;i, . . . , Zp\a) 
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= 5Z n{c"'^(A^-^.M|S;+(m-l)-«)y;(a+(m-l)-d)} f[ j^^" 

^01 7,1 "l^crfm) ■^a-(l) ) 

crGSp m=l l=m+l ^ ^ ' 

(50) 

We are now in the position to proceed to the general sl{n) case. 

We start with a special term of (0), where all operators T„_iq, are to the left of operators T„_i^ 
if a < jS. The associated coefficient $ again remains invariant under the action of F, which 
enables us to express it in a manner similar to (^^ . We then use the Faddeev-Zamolodchikov 
algebra to reverse the order of all operators. Once again only the term containing one ^ in 
the expression (^) of the respective operators Tn-ia contributes in this special ordering. 
The wavevector can be expressed in a form similar to (^) {a = n — 2 — a) 

n-2 Pa 

rpii-]) r 

iiT^.-.^jpg a=0j=pa+i+l 

(51) 



^„(iv,po,pi,...,Pn-2)= <.m,...,p.-.(^i'---'V+....„-.kn,---,^«.oia)n n ^?ii^5"^ 



with the recursion relation for B^"^^ 

-^POPI •••Pn-2('^1' ■ ■ ■ ' -^P0+Pl + ---Pn-2 • • • 5 ^Po\^) 

n-3 n—2 Pa P/3 

= EH n n n ka.(w-..j 

(o-l 
K(lf,) - K{ka)\0'+ J^(Pm -Pm+l) ■ m + [K-Pa+l - 1) " « + {1(3 - Pp+l " 1) " 
m=0 

n-2 n-2 

X n ^p^'-pI+i (Aa(p,+i+i) • • • A<,(p^) . . . 2;p^ |a + ^ (p„-, - p^+i) ■ m) 

7=0 m=7+l 

X -^Pi...pL2('^po+i • • • Apo+pi+...+p„-2|Act(i) • • • A<j(pi)|a) (52) 
which can be solved explicitely to yield 

-^POPI ■■■Pn-2 ("^1' ■ ■ ■ ' Apo+pi+...+Pn-2 1^1' • • • ? ^PO 1*^) 

n— 3 n—3—i n—2—i Pa^ P0^ 

(ToGSpp (TiG5p-^ (t„_3GS'p^_3 i=0 ai=0 j3i=ai+l fcaj=Pa^+i+l /^^=p^^^j^+l 

/ Oii-l 

I A,,_i+^,(,^j - Aq,_i+^,(fc,j|a + ^ (pm - Pm+i) ■ rh + (A;„^ - - 1) ■ di + (/^ - p^- +i - 1) ■ A 

\ m=0 
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n—2—i 

Ag,_2+a,_i(p^-j|a + ^ {Pm- Prn+l) " m) 

(53) 

where we defined 

i 

j=0 

and 

By expressing the sl{n) wavevector (|5l| ) with the help of sl{2) building blocks (|53|) we have 
achieved a resolution of the hierarchy. 

Remark : Once again the arguments in this section apply to the sl{n) XXZ-model as well, 
leading to a wavevector for this model which displays the same features as ( pTD in connection 
with (0), without the a dependence and with a trigonometric parametrization. 



n—2—% 



7«=0 



^ (Ag„-3 + l • • • Ku- 



3+Pn-2 I A(J„_4+o-„_3(l) • • • ^g„_4+(T„_3(p„_2) 



...A„ 







5 Conclusion 



We accomplished the construction of a factorizing F-matrix for the Belavin model enabling one 
to construct completely symmetric creation operators which moreover are devoid of non-local 
effects from polarization clouds. These operators were used to resolve the intricacies of the 
nested structure of Bethe wavevectors. 

In contrast to the sl{n) XXX-model the above method does not rely on an sl{n) invariance 
of the monodromy matrix, which renders the extraction of the elements in the grid of the 
monodromy matrix much more cumbersome. The only ingredients needed in our computation 
are the form of the R-matrix (^, i.e. its structure, unitarity and the fact that it constitutes 
a representation of the permutation group, and the property (|^). Thus our findings directly 
yield the corresponding expression for the generalized sl{n) XXZ-model. 
In view of the similarities between the generalized sl{n) rational, trigonometric and elliptic 
model it is tempting to ask whether such a distinguished basis exists for all integrable models 



in two dimensions. In |19| it was shown that in every integrable two-dimensional quantum field 
theory there exist semi-local polarization-free generators which are localized in wedge-shaped 
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regions of Minkowski space. It is conceivable that there is a relation between these operators 
and polarization-free operators in lattice spin models. 

We hope that these results might prove useful for the construction of formfactors starting 
from the microscopical level. 



Acknowledgement: We thank H. Boos, R. Flume and R.H. Poghossian for helpful discus- 
sions. 
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